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I : Abstract 

m 
(N 

^ . The linear response of a thin superconducting strip subjected to an appHed 

o . 

perpendicular time-dependent magnetic field is treated analytically using the 
, method of matched asymptotic expansions. The calculation of the induced 

0> , current density is divided into two parts: an "outer" problem, in the middle 



of the strip, which can be solved using conformal mapping; and an "inner" 
problem near each of the two edges, which can be solved using the Wiener- 



O . Hopf method. The inner and outer solutions are matched together to produce 

CJ 

a solution which is uniformly valid across the entire strip, in the limit that 
the effective screening length Aeff is small compared to the strip width 2a. 
From the current density it is shown that the perpendicular component of 
the magnetic field inside the strip has a weak logarithmic singularity at the 
edges of the strip. The linear Ohmic response, which would be realized in 
a type-II superconductor in the flux-flow regime, is calculated for both a 
sudden jump in the magnetic field and for an ac magnetic field. After a jump 
in the field the current propagates in from the edges at a constant velocity 
V = 0.772D/d (with D the diffusion constant and d the film thickness), rather 
than diffusively, as it would for a thick sample. The ac current density and the 
high frequency ac magnetization are also calculated. The long time relaxation 
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of the current density after a jump in the field is found to decay exponentially 

with a time constant tq = 0.255ad/D. The method is extended to treat 

the response of thin superconducting disks, and thin strips with an applied 

current. There is generally excellent agreement between the results of the 

asymptotic analysis and the recent numerical calculations by E. H. Brandt 

[Phys. Rev. B 49, 9024; 50, 4034 (1994)]. 
PACS numbers: 74.60.-w, 74.25.Nf, 02.30.Mv 
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I. INTRODUCTION 



When a thin superconductor is placed in a perpendicular magnetic field there are large 
demagnetizing fields which produce an enhanced response to the applied field. For instance, 
the induced magnetic moment in the perpendicular geometry for a sample with thickness d 
and a width a ^ c? is 0{a/d), while in for a longitudinal magnetic field it is only 0(1) [|T|. 
These demagnetizing fields are important not only for static properties, but also for dynamic 
properties, such as the response of the sample to an applied time-dependent magnetic field. 
Surprisingly, there has been little theoretical work on this subject until recently. Brandt 
has devised an efficient numerical method for calculating the linear or nonlinear response of 
superconducting strips or disks in time-dependent magnetic fields. His studies have 
unveiled a number of interesting properties of the sheet current and magnetic moment in 
these geometries, including dynamic scaling properties of the current density at the sample 
edges at high frequencies or short times. 

The present paper is an analytic treatment of the linear electrodynamics of thin supercon- 
ducting strips and disks in perpendicular magnetic fields, and as such is complementary to 
Brandt's numerical work. Many of the features of the linear response which were extracted 
numerically by Brandt emerge naturally from this analysis; the analytic results obtained 
here agree quite well with Brandt's numerical results. The calculational technique employed 
here is the method of matched asymptotic expansions a technique originally devised 

to treat boundary layer problems in fiuid mechanics f^, and used recently to study several 
interesting problems in nonequilibrium and inhomogeneous superconductivity P,^. The 
idea is to split the problem into two pieces; an "outer problem" in the middle of the strip 
(which is straightforward to solve), and an "inner problem" near each of the two edges of 
the strip (which requires a bit more ingenuity). The solutions are then matched together, 
and a "uniform" solution, valid across the entire strip, is constructed. The time dependence 
of the applied field is treated using Laplace transforms, which provides a unified framework 
for treating the transient response after the applied field is suddenly switched on, as well as 
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the steady-state ac response. It should be noted that there is an allusion to such a matching 
procedure for the same problem in a paper by Larkin and Ovchinnikov [|^; these authors 
simply quote a result for the behavior of the current near the edge of a strip in a perpendic- 
ular field. The present work goes well beyond that of Larkin and Ovchinnikov, by explicitly 
constructing the inner, outer, and uniform solutions, and using these solutions to study the 
nonequilihrium response. It also appears that the result quoted by Larkin and Ovchinnikov 
is incorrect in detail (see Appendix B). 

As this paper is somewhat long, the primary results are collected in Table |l[ The organi- 
zation of the paper is as follows. In Sec. II the integro-differential equation for the current 
density in the strip is derived, and a small parameter e = 2Aefr/a is identified, with Agfr the 
effective penetration depth and 2a the strip width. Sec. Ill treats the large-e limit, which is 
helpful for anticipating some of the features of the small-e solution. The asymptotic analysis 
for small-e is constructed in Sec. IV; the outer problem is essentially solved by conformal 
mapping, while the inner problem is solved using the Wiener-Hopf method. The solutions 
are matched using formal asymptotic matching, and the uniform solution is constructed. 
The uniform solution is used to study the perpendicular magnetic field in Sec. V, and the 
current density and magnetization in Sec. VI. In Sec. VII the method is extended to treat 
the response of thin superconducting disks, and thin strips in the presence of an applied 
current. In Sec. VIII the long-time behavior is studied via an eigenfunction expansion of 
the current density, and the fundamental relaxation time is calculated using a variation on 
the analysis of the previous sections. The results are summarized in Sec. IX. Some of the 
more algebra intensive parts of the calculation are relegated to Appendices A-C. 

II. DERIVATION OF THE INTEGRO-DIFFERENTIAL EQUATION 

To begin, we will derive the integro-differential equation which determines the current 
distribution in the strip. The geometry is illustrated in Fig. |I]. The strip is in the x — y 
plane, being infinite in the ^/-direction, and having a width 2a such that the strip occupies 
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the region —a < x < a. The apphed field Ha(t) = Ha{t)z is normal to the strip and in the 
z-direction. The vector potential A satisfies 

- V^A = 47rJ, (2.1) 

in the transverse gauge in which V ■ A = 0. We will focus here on situations in which the 
current density is invariant along the y-direction (along the length of the strip), so that 
both J and A are along the ^/-direction. In the thin film approximation, the current density 
Jy{x,z,t) is averaged over the thickness d of the film; the averaged current will be denoted 

by j{x, t), so we have 

Jy{x,z,t) = dj{x,t)S{z)e{a^ - x^). (2.2) 

Eq. ( p^.l| ) is solved by introducing the Green's function for the two dimensional Laplacian, 

G{x — x', z — z'): 



Ay{x, z,t) = Ao^y — An J G{x — x' , z — z')Jy{x, z,t) dx' dz' 

= Ao^y - And r G{x - x' , z)j{x', t) dx' . (2.3) 

J ~a 

Differentiating both sides with respect to x, and using dAQ^yjdx = Ha{t), we obtain 

dAy{x, z, t) 



H^{x,z,t) 



dx 

r« dG{x — x', z 



Ha{t) - And / ^jOr , t) dx' . 2.4 

J -a OX 



Finally, if we specialize Eq. ( p.4|) to 2; = 0, and use dG{x — x',0)/9x = \/2n{x — x'), we 
obtain for the magnetic field normal to the strip, 

H,{x,t) = Ha{t) + 2d{P) r ^-^^dx', (2.5) 

J —a Jy dj 

where (P) indicates a principle value integral. 

To complete the description, we require a constitutive relation between the averaged 
current and the fields. This paper will concentrate on the linear response of the current, so 
that the general time-dependent response is 



j{x,t)=f a{t-t')Ey{x,t')dt', (2.6) 

with a the conductivity. In most of this paper we shall be interested in the solution of the 
initial value problem; i.e., the time evolution of the sheet current after an applied current 
has been switched on. Then it is natural to Laplace transform the currents and the fields 
with respect to time t: 



oo 

-St 



Jix,s)= / e-''j{x,t)dt, (2.7) 







and so on. The inverse Laplace transform is 



jM = —I e''j{x,s)ds, (2.8) 

Zm Jc-ioo 



with the integration contour chosen to pass to the right of any singularities. After Laplace 
transforming, Eq. ( |2.1CI| ) becomes 

j{x,s) = a{s)Ey{x,s) 

= -sa{s)Ay{x,s), (2.9) 

where Ey{x, t) = —dAy{x, t) / dt and an integration by parts has been used to obtain the last 
line. Now Laplace transform Eq. (|2.5|) with respect to time, and use Eq. (|2.9| ) to eliminate 
the fields in favor of the currents, to obtain the following equation of motion for the current: 

-AnKs(s)d^^^ = Ha{s) + 2diP) r ^-^y^dx', (2.10) 

ox J -a X' — X 

where Aefr(s) is an effective screening length defined by 

For a superconductor, (t(s) = l/47rA^s, with A the London penetration depth, so that in this 
case Aeff = A^/rf. For an Ohmic conductor, Acfr = l/47rcr(0)c?s = D/ds^ with D the diffusion 
constant for the magnetic flux. An Ohmic response would be realized in a normal metal or 
in a type-II superconductor in the fiux-fiow regime. Eq. (|2.10|) is essentially identical to the 
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equations of motion derived by Larkin and Ovchinnikov Eq. (33), and Brandt |T[], Eq. 
(3.6) the only difference being that Laplace rather than Fourier transforms are being used 
here. This equation of motion can be conveniently expressed in dimensionless variables by 
writing x' = x/a, f{x', s) = j{x, s)/{Ha{s)/ 271(1), and e = 2Acfr(s)/a: 

- efix, s) = l + -(P) ^^^^ dx' (2.12) 

TT J-1 X' — X 

(the s dependence of e will be suppressed for notational simplicity; it will be reinstated later 
when the transforms are inverted). The current density is the primary quantity of interest 
in this paper. A related, and experimentally accessible, quantity is the magnetic moment 
per unit length of the strip 



M{s) = I xj{x,s)dx 

—a 



Ad 

where mis) is a dimensionless moment defined as 



m{s) (2.13) 



2h{s) 

m(s) = — xf{x,s)dx, (2.14) 



7T 



where h{s) is defined through Ha{s) = Hah{s). 

There are no known analytical solutions of Eq. ( p.l2| ); its solution is the subject of the 
remainder of this paper. However, for a typical sample e <C 1, in which case the left hand 
side of Eq. ( p.l2| ) constitutes a singular perturbation (the highest derivative is multiplied by 
the small parameter), and we can bring to bear all of the techniques of asymptotic analysis 
to solve this problem . Before attempting this, we will first develop a perturbative analysis 
which is valid for e ^ 1, which is physically less interesting but mathematically simpler. 



III. EXPANSION FOR LARGE e 



To study the behavior of Eq. ( p.l2| ) for large e, it is useful to first rescale by introducing 



a new function g{x, s) = ef{x, s). We then expand g{x, s) in powers of e ^: 



g{x,s) go{x,s) + e gi{x,s) + .... 



(3.1) 



Substituting into Eq. ( p. 12 ) and matching terms of the same order, we have 



9oix, s) = -1, 



i(p) r'^-^dx' 

TT 



(3.2) 



(3.3) 



1 X' — X 

and so on with the higher order terms. Assuming that there is no net current in the strip 
(i.e., no apphed current), then the current is odd in x, so that g{0) = 0, and we have 



go[x,s) = -X, 



gi[x,s) 



1 

7r 



1 — , fl + x 
X H : — m 



(3.4) 



(3.5) 



2 Vl-x. 

The current increases hnearly across the sample, except near the edges. Near the left edge 
{x = —1), we have 

1 



g{x,s) ~ 1 



7re 



[l + {l + x) ln(2/(l + a;))] . 



(3.6) 



As we shall see below, a similar behavior near the edge will also emerge in the small-e limit. 
From the expansion for f{x, s) we can calculate the magnetic moment from Eq. ( |2.14| ): 

"4 1 2 



mis) 



his) 



+ 0{e- 



(3.7) 



3 7re(s) (7re(s))^ 

For an Ohmic conductor, e(s) = 2D /ads, so the magnetic moment vanishes linearly with 
frequency at low frequency. This result is essentially equivalent to Eq. (5.6) of Ref. JIl. 



IV. ASYMPTOTIC ANALYSIS FOR SMALL e 



We now turn to the solution of Eq. ( p.l2[ ) for small e. An asymptotic solution can be 



obtained using the method of matched asymptotic expansions. We break up the strip into 
an "outer region," which is the interior of the strip, and two "inner regions," one near each 
edge. This is illustrated schematically in Fig. 0. The solutions are then matched in common 
overlap regions. A detailed discussion of the method can be found in Ref. 0]. 
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A. Outer solution 



The expansion in the outer region (away from the edges) is obtained by expanding / as 
a series in e: 

/(x;e)~/o(x) + e/i(x) + ... . (4.1) 

Substituting this expansion into Eq. ( p.l2| ) and collecting terms of the same order, we find 
for /o(a;) 

\p) f' M^dx' = -1. (4.2) 



The solution of this singular integral equation which is odd is x can be found in Ref. 0: 

/o(^) = - (i_^2)i/2 ' (4-3) 

which coincides with the usual solution obtained from conformal mapping techniques. Near 

the edges at x = ±1, /o behaves as 

SO the current has a square root divergence at the edges. This is due to the fact that the 
outer solution corresponds to complete screening of the applied field, which can only be 
achieved by having an infinite current density at the edges. The outer solution therefore 
breaks down at distances of order e of the edges; the current at the edges is thus of order 
e^^/^. In order to remedy this problem, we proceed to the solution of the inner problem near 
each of the edges. 

B. Inner solution 

We will first study the inner problem at the left edge, x = —1. The outer solution breaks 
down at X + 1 ~ e, suggesting that the appropriate variable in the inner region is 

X = {x+l)/e. (4.5) 



Also, since /(—I) ~ e we rescale f(x) in the inner region as 

F{X)=e'/'f{x), (4.6) 
so that -F(O) = 0(1). In terms of these inner variables, Eq. ( p.l2[ ) becomes 



F' = -e'/^-l{P)l'''^gldX'. (4.7) 

Next, expand F{X;e) in powers of e^^^, as suggested by the rescaled form of the integro- 
differential equation: 

F{X- e) ~ Fo(X) + e^/2Fi(X) + . . . . (4.8) 
The lowest order term satisfies 

For small e the second integral will be dominated by the large X behavior of Fq{X); it will 
be shown below that in order to match onto the outer solution this is necessarily of the form 
Fo(X) ~ (2X)-i/^ Therefore, we see that the second integral is of order e^/^, and can be 
dropped at this order of the calculation. Our final integral equation for -Fo(-^) is then 

The problem in the inner region consists of solving a homogeneous integro-differential equa- 
tion on a semi-infinite interval. This is equivalent to finding the current distribution in a 
semi-infinite strip in zero applied magnetic field. 

Eq. ( 4.1U| ) can be solved using the Wiener- Hopf method as follows. The function 



Fo{X) — s> as X ^ oo, and Fo{X) = for X < 0; introduce a second unknown function 
G'(X) such that G{X) = for X > and G{X) ^ as X ^ -oo. We then introduce the 
complex Fourier transforms of -Fo(X) and G{X), 

roo 

$+(A;) = / Fo(X)e^'=^dX, (4.11) 
Jo 
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G-ik) = f G{X)e'^^dX, (4.12) 

J —oo 

such that is analytic for Im(A;) > — /? and G^{k) is analytic for Im(A;) < a, for some 

a > (3. Then Fourier transforming Eq. ( |4.10| ), and integrating by parts, we obtain 

G_(A;) - Fo(0) - ik^+{k) = isgn{k)^+{k). (4.13) 

We see that as /c — >• oo, 

<^>+(k) -^^ + 0{k-^). (4.14) 
ik 

To take care of the ambiguities in defining sgn(fc), replace it by /c/ (/c^ + 5^)^^^, with the real 
part > for Iie{k) > 0, and choose the branch cuts to run between (—zoo, —i6) and {i6, too). 
We can then take 5 — at some later point in the calculation. Rearranging Eq. ( [4.13| ) a 
bit, we have 

ikK{k)^+{k) = G-{k) - Fo(0), (4.15) 

where 

Now, if we can factor K{k) into the form K{k) = K^{k)/ K^{k), with K+{k) analytic for 
Im(/c) > —6 and K_{k) analytic for Im(/c) < 5, then we may rewrite Eq. ( [4. 151) as 

ikK+{k)^+{k) = K_{k) [G-{k) - Fo(0)] . (4.17) 

Both sides are now analytic in their respective regions of analyticity; we can then use analytic 
continuation arguments to note that both sides must then equal an entire function E{k). By 
examining the limiting behavior of the left hand side as A; ^ oo, we see that this function 
must be chosen to be a constant G (any positive power of k would produce non-integrable 
singularities in Fq{X)), so that we finally have 
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The function Fo{X) is then obtained by inverting the Fourier transform, 



oo 



Fo(X) = ^ I ^+{k)e-''''dk, (4.19) 

with the integration path indented so as to pass above any singularities on the real axis. 

The only remaining task is the decomposition of K{k), which is carried out in Appendix 
A. Using Eqs. ( |A^ ) and ([ASD, we have 

e '"^ dk 



OO A;(l + l/|A;|)i/2 
C ['^ ?,m[kX + ^{k)] 



7r Jo (F + fc)V2 

where the phase ^{k) is given by 



dk, (4.20) 



m(k) = - + - -du. (4.21) 

The constant C is determined from the matching conditions, which are discussed below. 



C. Asymptotic matching 

The inner and outer solutions may now be matched together in a suitable overlap region. 
This is done by expressing the outer solution fo{x) in terms of the inner variable X = 
{x + l)/e, and then taking X —>■ while holding e fixed: 

eX - 1 
[(2-eX)eX]V2 

^ (X ^ 0). (4.22) 



fo{X) 



(2eX)i/2 

The inner solution Fo{X) must match onto this outer solution as X ^ oo, so the asymptotic 
behavior of -Fo(-^) niust be 

i^o(X)~^^^ (^-oo). (4.23) 

Now expand Eq. ( [4.2(J| ) for large-X; the integral is dominated by the small-fc behavior of the 
integrand, and we find 
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SO that we have 

C = - (4.25) 

in order to match the inner and outer solutions. Therefore our final expression for the inner 
solution is 

p^(X)- ^ f^ sm[kX + ^{k)] 



Comparing Eqs. and ( jO^ ), we see that Fo(0) = -C = (7r/2)i/2. 



By rotating the integration contour, it is possible to show that for X < the integral 
vanishes (as it should), while for X > the integral may be written as 

\ poo p-Xy-g{y) 

^»w = (2^i + 1)3/. 

with 



9{y) = -r^,du. (4.28) 

TT JO 1 + 



This function has the limiting behaviors 

-2/ln(e/y)/7r + 0(y3), y < 1; 
9{y) = { (4.29) 
-ln(ei/)/7ry + 0(y-3), y > 1. 

The square-root singularity in the integrand at ?/ = can be removed by changing variables 
to y = z^, so that 

.9x1/2 .00 ^-Xz'-g{z^) 

which is particularly convenient for numerical evaluation. As shown in Appendix B, for 
small-X, -Fo(-^) has the expansion 

Fo{X) = —^rrj^ In -1A228X + XlnX + 0{X^)] (4.31) 



(27r)i/2 
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(which disagrees with Eq. (37) of Larkin and Ovchinnikov, Ref. 0; see the remarks in 
Appendix B). Although the current is finite at X = 0, it has a slope which diverges as 
ln(X). The result of a numerical evaluation of the integral, along with a comparison of the 
numerical results to the asymptotic expansions, is shown in Fig. ^ 

So far we have only discussed the matching procedure at the left edge of the strip x = —1. 
The same procedure can be carried out at the right edge, x = 1, as follows. At the right 
edge, the inner variable will be X = (1 — x)/e. As before, we will also rescale f{x) as 
F{X) = t^/^f{x). Substituting these expressions into Eq. ( p.l2| ), we obtain 

F\X) = e'l' - i(P) fj' 1^ dX\ (4.32) 

which is the same as Eq. ( |4.7|) except for the minus sign in front of the e^^^. Expanding 
in powers of e^/^, the 0(1) term, Fo(X), satisfies Eq. (|2.12|) , and the method of solution is 
identical. To match onto the outer solution, we write fo{x) in terms of X, and then take 
X ^ while holding e fixed: 

1-eX 



fo{X) 



[(2 - eX)eX]V2 

^ (X ^ 0). (4.33) 



(2eX)V2 

Therefore, the asymptotic behavior of Fq{X) must be 

Fo{X) ~ -^^^ {X ^ oo). (4.34) 

Therefore, Fq{X) = —Fq{X), which could have been surmised from the symmetry of the 
problem. 



D. Uniform solution 

We are now in a position to construct an asymptotic solution which is uniformly valid 
across the entire width of the strip; i.e., valid for all x as e — 0. To do this we simply add 
the inner and outer solutions; however, this would produce a result which was 2/match(a;) in 
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the matching region, so we also need to subtract the /match (a^) for each of the two matching 
regions [Q. The result is 

f ( \ - ^ M ^ ^ 

Junif l^a;, Sj 9M /o ' 



(l-x2)l/2 [2(l-x)]l/2 [2(1+X)]V2 

+e{s)-'/' {Fo[(l + x)/eis)] - Fo[(l - x)/e{s)]} . (4.35) 

The uniform solution is plotted in Fig. ^ for e = 0.1. 

With the uniform solution we can calculate the magnetic moment, given in Eq. ( p. 141) . 
The result is 



m 



^ hf^^^ 8 , f2\'/' ^oo e-g(-^)[(y + l)e-^^ + y - 1] ^ ] 



. 2 ^^ 3/2 ^oo e-giey) [(^^ ^ l)g-2j; + y _ l] 
|/5/2[l + e2y2]3/4 

For e = 0, the integral is 8/37r, so that m(s) = h{s), which is the ideal screening limit. 
Determining the leading e behavior of the integral is rather subtle; the details are relegated 
to Appendix C. The result is 



m(s) = h{s) 



(4.37) 



This expression should give the correct high-frequency (s oo) behavior of the magnetiza- 
tion. Note that for s = (e ^ oo), m{s) = 1 — 8/37r = 0.151; however, we know that the 
correct limiting behavior is m(0) = (see Sec. III). Therefore the uniform approximation 
does not reproduce the correct low frequency behavior of the magnetization for an Ohmic 
conductor. 



V. MAGNETIC FIELD WITHIN THE STRIP 



By using the constitutive relation, Eq. (|2.9|) , it is also possible to calculate the magnetic 
field perpendicular to the strip. Using Hz{x, s) = dAy{x, s)/dx, going to our dimensionless 
variables, and using the uniform approximation from the section above, we have 

df{x,s) 



Hz{x, s) 



Ha(s) le(s) 



dx 



1 



1 



-e s 



-1/2 



(l-a;2)3/2 [2(l-a;)]3/2 [2(l + x)]3/2 
{F^[(l + x)/e(s)]+F^[(l-x)/e(s)]}}. 
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(5.1) 



The magnetic field is plotted in Fig. ^ From the results in Appendix B, close to the edges 
(1 ± X = 0(e)) we have 

KiX) = HyX) + 0(1), (5.2) 

so that it would appear that the field diverges logarithmically at the edges, in agreement with 
the numerical work in Ref. However, as we get even closer to the edge (l±a; = 0(e^)), this 
log divergence is swamped by a square-root divergence from the outer and overlap terms. 
This latter behavior is most likely an artifact of the approximation, and would probably 
disappear in a higher-order calculation. 



VI. DYNAMICS OF THE CURRENT DENSITY AND MAGNETIZATION 

Having obtained a uniformly asymptotic solution to the equation of motion for the av- 
eraged current density, we can now examine its evolution in the time domain by inverting 
the Laplace transform for j{x, s), Eq. ( p.8| ). The details of the inversion process will depend 
upon the time dependence of the applied field, and the model chosen for e(s). Two different 
models for e(s) will be considered: (1) e(s) = /ad a constant, corresponding to a super- 
conductor; (2) e(s) = {2D/ad){l/s), with D = l/47rcr(0) the diffusion constant for flux in 
the normal phase, corresponding to an Ohmic conductor (a type-II superconductor in the 
flux- flow regime, for instance |jl|). 



A. Superconductor 



For a superconductor the inversion of the Laplace transform is particularly simple. In 
conventional units we have (recall that Aefr = X^/d) 

Ha{t) 



Junif i^yt) 



2nd 



X 


1 


1/2 


1 


1/2 


(l-x2)l/2 + 


2(1 -x)_ 




_2(l + x)_ 





'1/2 



1 + x 



1 — a; 



e / V e 

In this case the induced current is in phase with the applied fleld. 
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(6.1) 



B. Ohmic conductor: penetration of a jump in the applied field 



We will first treat the case in which the field is suddenly switched on, so that Ha{t) 
Ha9(t), and thus Ha{s) = Ha/s. Inverting the Laplace transform, we have 



Junif t) 



Ha{t) 



2nd 



(1 -x2)l/2 



nl/2 



2 1 -X 



2 1+X 



1/2- 



+jedge(l + X,t)- jedge(l -X,t), 



(6.2) 



where the edge current ieAge{x,t) is given by 



1/2 



jedge (-^j ^) 



dy 



Hg ( ad \ 

2nd\W) Jo yl/2(^/2 + 1)3/4 

1 /-c+ioo Qxp[st - {adx/2D)ys] 



ds. 



(6.3) 



Zm Jc—ioo 

(7rs)i/2 

The Laplace transform can be calculated by closing the contour in the left half plane (for 
t > dx/2D), wrapping the contour around the branch cut along the negative s axis, with 
the result 



Ha f ad ^ / '^'^ ^\ 

Jedge(x,t) = ^^ :F^i—- , 



(6.4) 



where J-'i{u) is a scaling function given by 



-9(2/) 



1 rV" 

^'^""^ ^nJo (1 - My) 1/2^1/2 (^2 + 1)3/4 



dy. 



(6.5) 



The scaling function has been normalized so that J^i(O) = 1; for large-w, J-'i{u) ~ m"^/^. 
For numerical purposes, it is useful to transform the integral by making the substitution 
1 — uy = sin^ 9, so that 

2u f^/"^ exp[—g{cos^ 6/u)] 



^i(w) 



de. 



(6.6) 



Vr Jo (cOs4^ + m2)3/4 

which removes the square-root singularities from the integrand at the endpoints of integra- 
tion. 

The scaling function J-'i{u) is plotted in Fig. ^ There is a maximum at u = 0.386, 
with JFi max = 1.1078. From the position of this peak we can define a velocity v of fiux 
penetration: 
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v=(-) = 0.772^. (6.7) 

\t / max d, 

These results agree exactly with the numerical work of Brandt As noted by Brandt, 

in the thin film geometry the flux entry is ballistic rather than diffusive, as it would be in a 
bulk sample with no demagnetizing fields. 

We can also calculate the time dependent magnetization after a jump in the field. By 
using the small-e expansion in Eq. (|4.371 ), and inverting the Laplace transform by wrapping 
the integration contour around the branch cut along the negative-s axis, we find at short 
times 

m{t) = 1 - —M/t) \n{8ne-^/H/r) + Oi^), (6.8) 

where r = ad/27rD is a characteristic relaxation time. Similar behavior was found by Brandt 
in his numerical studies of the magnetization; for the prefactor of the log he obtained 
0.205, compared to the present value of G/tt^ = 0.194; for the constant inside the log, he 
obtained 25, compared to our value of 8776"^/^ = 12.9. It is not clear whether these small 
discrepancies are the result of the approximations in this paper or uncertainties in Brandt's 
numerical work. 



C. Ohmic conductor: ac response 

Next, we consider the response of the strip to an ac magnetic field, Ha{t) = Haexp{iujt), 
so that Ha{s) = Ha/{s — iuj). When inverting the Laplace transform, there will be contri- 
butions both from the pole at iuj and from the square-root branch cut along the negative-s 
axis. The branch cut contribution decays as t~^/^; since we are concerned here with the 
steady-state behavior, we will neglect this term, keeping only the pole contribution. Writing 
jedge{x,t) = jedge{x,u) exp(zco't), wc havc for the current near the edge 

Ha f LJ7Tad\^^^ ( dujx\ 

Jed,e(x,^) = — (^^j (6.9) 

where the ac scaling function J-'2{u) is given by 
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-g{y)-iuy+i-K/4: 



■dy. 



(6.10) 



TT Jo ?/l/2(y2 + 1)3/4 

The scaling function is defined so that ReJF2(0) = ImjF2(0) = 1. The real and imaginary 
parts of the ac scaling function are plotted in Fig. |^. The real part of the scaling function 
has a maximum at m = 0.232 {xuo = 0.0738 a/r in conventional units), with RejF2iiiax = 
1.0787; the imaginary part changes sign at m = 2.43 {xu = 0.773 a/r in conventional units). 
These results are once again very close to Brandt's numerical results The uniform 

approximation to the current density is 



2iTd 



X 



(1 _ 3-2)1/2 



+ 



2(1-0;) 



1/2 



1 



2(1 + x) 



1/2^ 



+jcdgc(l + X,U) - Jedge(l - X,u). 



(6.11) 



Using the small-e expansion of the magnetization, Eq. ( [4.37| ), we can also calculate the 
high frequency magnetization; the result is 

6 ln(87re^~^/32cjr) 



m[u] 



lUOT 



-2\ 



(6.12) 



which is quite similar to the numerical result obtained by Brandt ]l|. The constants differ 
slightly; here the prefactor is 0.194, while Brandt obtained 2/7r^ = 0.203; for the constant 
inside the log, Brandt obtained 16.2, compared to the present value of ^Tie^'^^'^ = 8.45. 



VII. EXTENSIONS OF THE METHOD 

With minor modifications it is also possible to treat two related problems, the current 
distribution in a thin superconducting disk in a perpendicular field, and the current distri- 
bution in a thin strip in the presence of an applied current. Rather than discussing these 
cases in detail, only a brief sketch of the results will be provided. 



A. Disk geometry 



Rather than a strip we now have a superconducting disk centered at the origin of the 
X — y plane, of radius a. The current density and the vector potential are both in the (f) 
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direction. Using a Green's function method similar to that in Sec. II, the 2;-component of 
the magnetic field inside the disk satisfies 

H,{r, z = 0, t) = Ha{t) + 2d{P) r P{r, r')j{r', t) dr', (7.1) 

Jo 

where the kernel P{r,r') is given by 

K(k) E(k) ,2 4rr' , , 

P{r,r') = ^ + ^, k'= (7.2) 

where K{k) and -E(/i;) are complete elliptic integrals of the first and second kind. We 
again Laplace transform, use the constitutive relation between j and A, and go to the 
dimensionless variables r' = r/a, /(r, s) = j{r, s)/ {Ha{s)/TT^d), and e(s) = 2Aefr/a, to arrive 
at the following integro-differential equation: 

-e-|-[r/(r,s)] = ^ + -(P) f P{r,r')f{r' , s)dr' . (7.3) 



r (9r 2 IT Jo 

The outer solution is 

which has a square root singularity at the edge. Near the edge we construct the inner 
solution by defining the inner variables R = (1 — r)/e, F{R) = e^/^/(r); for small e the 
kernel behaves as 

P(l -eR,l- eR') = -^^^ + O(lne). (7.5) 

After expanding F{R] e) in e, we find that the lowest order term satisfies Eq. (|4.1CI| ), so that 
the inner solution is the same as before. Finally, the inner and outer solutions are matched 
as before; the uniform solution is then 

/-f(r, s) = ^^37^ - [2(1 - r)]i/^ + <^r'^'Fo[il - r)/e{s)]. (7.6) 
The magnetic moment of the disk is 

/•a 

M(s) = 71 r'^j(r,s)dr 
Jo 
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where the dimensionless moment in the disk geometry is 



m{s) = r^f{r, s) dr. (7.8) 



2 Jo 

Using the uniform approximation, this becomes 

m(s) = h(s) < 1 = / ^ > . (7.9) 

For small e the integral can be expanded using the same method as for the strip (see Ap- 
pendix C), with the result 



m(s) = h{s) 



(7.10) 



For an Ohmic conductor e(s) = 2D/ads = l/nrs, with r = ad/2'KD. After inverting the 
Laplace transform, we find that after a jump in the magnetic field, the magnetization for 
small if: is 

mit) = 1 - ^(t/r) In Une-^^'H/T) , (7.11) 

while the ac magnetization at high frequencies is 

, , ^ 2v/21n(47re^-5/3^cjr) 
m{uj) = l-—^ ^ : -. (7.12) 

The behavior is essentially the same as for the strip, with slightly different constants. Similar 
behavior has been found by Brandt p| in his numerical studies. He finds a prefactor of 
S/tt^ = 0.304, compared to the present value of 0.286; for the constant inside the log, he 
obtains 11.3, compared to our 4.23. As in the strip geometry, the source of the discrepancy 
is unclear. 

B. Current distribution in the presence of an applied current 

The effect of an applied transport current is to modify the 0(1) outer solution, which 
now becomes 
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where fa = I{s)/{Ha{s)a/2), with /(s) the total current (not the current density) in the 
strip. Writing the outer solution in terms of the inner variable X, we see that the matching 
conditions are 

/o(X) = ±(l±/J^^-^, (7.14) 

with the + corresponding to the left edge and the — to the right edge. The solution to the 
inner problem is the same as before. After matching the inner and outer solutions, we have 
for the uniform solution 

X- fa , 1 - /a 1 + /a 



/'unif (2^) s) 



(l-a;2)i/2 [2(l-a;)]i/2 [2(l + x)]V2 

+<sr'^' {(1 + fa)Fo[{x + l)/e{s)] - (1 - /jFo[(l - x)/e{s)]} . (7.15) 



With this expansion it is possible to study the time-dependent response, just as in the zero 
current case. 



VIII. LONG TIME BEHAVIOR AFTER A JUMP IN THE APPLIED FIELD 

All of the previous sections have been concerned with the short-time or high frequency 
response of a strip to an applied field. In this section we will treat the long-time relaxation 
of the current density in an Ohmic strip after a jump in the perpendicular field. We will 
again use the method of matched asymptotic expansions, but in a slightly different form. 

We start with the equation of motion for the magnetic field, Eq. ( p.5|) , and differentiate 
with respect to time t. The time derivative of the magnetic field can be related to the 
current density for an Ohmic conductor with conductivity cr through dH^/dt = —dEy/dx = 
—a^^dj/dx. Since the applied field is a step function, dHa(t)/dt = Ha6{t). Therefore, for 
t > the current density satisfies (using a as the unit of length) 

- = 2rf(p) r dx'. (8.1) 

aa ox J-i x' — X ot 
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Following Brandt write the current density as an eigenfunction expansion: 



j{x,t) = ^^c„^„(x)e-*/^", (8.2) 



where the the relaxation times r„ = ad/2DXn are related to the eigenvalues A„, which follow 
from the solution of 



dx 



^(P) f'tA^dx'. (8.3) 

71 7-1 x' — X 



This equation is similar to the homogeneous version of the integral equation for the current 
density, Eq. (|2.12|) , but with an important sign difference on the left hand side. As a result, 



we can expect the eigenf unctions to be oscillatory, rather than decaying, in the middle of 
the strip. The long time behavior of the current density will be controlled by the smallest 
eigenvalue Aq, which produces the longest relaxation time Tq. 

We can develop an asymptotic analysis of the eigenvalue spectrum by first assuming that 
A„ ^ 1, so that 1/A„ serves as our small parameter. The consistency of this assumption 
should be checked at the end of the calculation. As before, we break the problem up into an 
outer problem in the middle of the strip, and two inner problems near each of the two edges. 
First we treat the outer problem. Define the outer variables = A^x, \l/^°''(Xo) = ipnix). 
Then the integral equation for the outer function is 

- -(^) / J" ° dX',. (8.4) 



dXo n J~\n X'^ — Xo 

Taking A„ ^ oo, we then have an integro-differential equation which relates the derivative 
of a function to its Hilbert transform. The solutions are cos(Xo) and sin(Xo); however, in 
the absence of an applied current the current density must be odd in x, so the physically 
acceptable solution is 

^i°)(Xo) = A„sin(X,), (8.5) 

with An a constant which can in principle depend upon n. This outer solution must be 
matched onto the inner solution, which we turn to next. 
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Let's first consider the inner problem at the left edge. Define the inner variables Xj = 
A„(H-a;) and "^^^Xi) = ipni^). Writing Eq. in terms of the inner variables, and taking 
A„ — cxD, we have 

^^-V)/ i^dX:. (8.6) 



dXi tt' 'Jo X[~Xi 

which is once again an integral equation of the Wiener-Hopf type. To solve, introduce an 

unknown function G{Xi) such that G{Xi) = for Xi > 0, and introduce the complex Fourier 
transforms 

poo 

= / ^^\X,)e'''''^dXi, (8.7) 
Jo 

G+{k) = / G{Xi)e'''^'dXi, (8.8) 



such that $+(fc) is analytic for Im(fc) > —(3 and Gj^{k) is analytic for Re(fc) < a, for some 
a > (3. Fourier transforming Eq. (p.6|), we then obtain 



G+{k) - v[/W(0) = ^kk{k) i^r^) K{k), (8.9) 



where 



U2 _ 1.2 



^.10) 



with 5 a small parameter which is taken to zero at some convenient point of the calculation, 
and kQ = [1 — 5^)^/^. The kernel K has been constructed so that it is free of zeros in the 
strip —5 < Im(A;) < 5, and K{k) ^ 1 as \k\ oo. The kernel can be factored into the 
quotient form K = Kj^/K^ using the general factorization procedure (see Appendix A), 
with the result that 

K+{k) = [K{k)]^/^e'^^''\ (8.11) 
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With the factorization, Eq. ( |8.9| ) can be written as 

{k - t5)k4k)[G4k) - ^«(0)] = ^k {j^^ K^m^{k). (8.13) 

Both sides are analytic in their respective regions of analyticity; analytic continuation allows 
us to set both sides equal to an entire function E{k). To choose E{k)^ we require that 
\l/^*^(0) be finite, so that ~ l/A; for large k. This can only be achieved by taking 

E{k) = Bnk/i^^"^, with Bn a constant (the 2^/^ has been added to simplify some of the 
resulting expressions). Therefore, 

Now we invert the Fourier transform, with the integration path passing above the poles on 
the real axis. Closing the contour in the lower half plane, we pick up contributions from the 
poles and from a branch cut which runs along the Im(A;) < axis, with the final result 

^t\X,) = -S„ cos(X, - 7r/8) + 5„^,,t(X,), (8.15) 

where ifjcntiXi) is the contribution from the branch cut, which is 0{X^^^'^) for large Xj. 
We must now match together our outer solution, Eq. (^.51), and our inner solution, Eq. 



g.l5| ). Taking the outer limit of the inner solution, and rewriting Xo and Xi in terms of x. 



we find that 

— Bn cos[A„(l + x) — 7r/8] = An sin(A„x). (8.16) 

This is only satisfied if A„ — 7r/8 = (n + l/2)7r and An = (— l)"'i?„. Therefore, the eigenvalue 
spectrum is 

A„ = — + n7r, n = 0,1,2,.... (8.17) 
o 

The same matching procedure must also be carried out at the right edge, and the analysis 
is identical. From the inner and outer solutions we can construct a uniform solution, which 
is 
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V^n,unif(x) = sin [{n + 5/8)7rx] + (-1)" {^cut[A„(l + x)] - V^cut[A„(l - x)]} , (8.18) 

where an overall constant has been dropped. 

This asymptotic expansion should be accurate for large n. To compare these results 
to Brandt's numerical work [llJ^, first note that Brandt calculates A„ = A„/7r, so we find 
Aq = 5/8 = 0.625. Brandt obtains the numerical value of Aq = 0.638, which is within 2% of 
the result of our asymptotic analysis. The approximation only improves for large n, so our 
result appears to be quite accurate for all n. This due in part to the fact that Aq = 57r/8 is 
large enough for the asymptotic analysis to be effective. With regard to the eigenf unctions, 
for large n the cut contributions become less significant (their contribution is localized near 
the edges), and so ipn{x) = sin [(n + 5/8)7ra:;] becomes an accurate approximation to the 
eigenf unctions. Based on his numerical results, Brandt quotes a similar result, but with 
5/8 replaced by 1/2; our result should provide a better approximation, and even appears 
to resemble the numerical result for n = 0. The cut contributions have derivatives which 
diverge logarithmically at the edges, similar to the numerical solutions. 



IX. DISCUSSION AND SUMMARY 

By applying the method of matched asymptotic expansions to the integro-differential 
equation for the current density in the strip, Eq. (|2.12| ), we have been able to derive a uniform 



approximation for the current density, which has been used to study the nonequilibrium 
response of the current in the strip, as well as the ac current density and the magnetization. 
Most of the effort has gone into understanding the response of an Ohmic strip. However, 
the method is easily generalized to more complicated dispersive conductivities ct(s); the only 
difficulties arise in inverting the Laplace transform to obtain the temporal response. For the 
purely Ohmic response, we found that after a jump in the perpendicular magnetic field the 
current propagates in from the edges at a constant velocity, in contrast to the longitudinal 
case, where the current propagates diffusively [|I| . The difference is due to the demagnetizing 
effects in the perpendicular geometry; initially the field lines must bend around the edges 
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of the sample, resulting in a large magnetic "pressure" which drives the current into the 
sample. 

A number of simplifying assumptions were made in order to make this problem analyti- 
cally tractable. The most important is the assumption of linear response. For many type-II 
superconductors, however, the current- volt age characteristics are highly nonlinear due to the 
collective pinning of the flux lines. In the perpendicular geometry there has recently been 
some progress in incorporating pinning (nonlinear response) into calculations of the current 
and field patterns for thin superconducting strips K is possible that the asymptotic 

methods used in this paper would be useful for studying the nonequilibrium, nonlinear re- 
sponse in the perpendicular geometry. A second assumption is that the current in the strip 
does not vary in the y-direction, so that we have an essentially one-dimensional problem. It 
would be interesting to include small variations of the current along the y-direction in order 
to determine the stability of the current fronts which enter after a jump in the perpendicular 
field; this problem might also be amenable to the type of analysis discussed in this paper. 
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APPENDIX A: DECOMPOSITION OF THE KERNEL K 

We want to decompose the kernel K{k) given in Eq. ( [4.16| ) into K^{k)/K^{k), with 
analytic in the upper half plane and analytic in the lower half plane. Since K{k) is free of 
zeros in the strip —6 < lm{k) < 6 and approaches 1 as ±oo. In K{k) is analytic in this 
strip and approaches as k ^ ±oo. We can therefore perform an additive decomposition of 

In K{k) = \nK+{k) - \nK_{k) (Al) 
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by writing 



In K+{k) 



1 roo-ia luK{z) 



2TTi 



z — k 



dz, 



(A2) 



with a > 5] there is an analogous expression for K^{k). This integral will be analytic for 
lm{k) > —6. Now if k can be taken to be real, and the integration path coincides with the 
real axis (indented to pass under the pole at k), then 

k f<^\n[K{x)/K{k)] 



1 k 
In K+{k) = -\nK{k) + — / 
Z m Jo 



x"^ — k^ 



(A3) 



where K{—x) = K{x) has been used to simplify the integral. Therefore, our decomposition 
ofK = K+/K^ is (letting 5^0) 

-1 ±1/2 



K^{k) 



(A4) 



with 



sgn{k) 



n Jo 



In 



u{\k\ + l) 



du 



v?-l 



(A5) 



\k\u + l 

The last integral may be rewritten in a more convenient form by first differentiating with 
respect to fc, integrating with respect to m, and finally integrating with respect to k: 



^(^) = — Z — /„ ^ — + TSgn(/i;). 



vr 



1 — 



(A6) 



This integral can be expressed in terms of dilogarithm functions, but this is not particularly 
useful for our purposes. 

For large |/c|, ^{k) has the asymptotic expansion 



1 ln(e ) 

TT k 

while for small /c, ({>{k) has the series expansion 



^{k) = ^sgn{k) - - ln(e/|A;|)A; + Oik^). 

4 TT 



(AT) 



(AJ 
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APPENDIX B: SMALL X BEHAVIOR OF Fo{X) 



In this Appendix the behavior of -fo(X) for small X will be derived by using the method 
of matched asymptotic expansions to derive a uniformly valid expansion for the integrand 
of Fo(X). Call the integrand in Eq. (g]^ I{y,X): 

^fa-^) = yV.(,. + 1)3/4 - (Bl) 

Expanding for small X, 

p-a{y) y^l'^p-9{y) 

This constitutes our inner expansion, which breaks down when Xy = 0(1). To derive the 
outer expansion, define an outer variable Y = Xy, rewrite I{y, X) in terms of Y, and expand 
to lowest order in X: 

^-Y~g{Y/X) ^ 
h(X, X) = yi/2^y2 _j_ X2)3/4"^ 

= ^X' + 0{X% (B3) 

In order to match the two expansions, express the inner expansion Ii{y,X) in terms of the 
outer variable Y, and expand for small X: 

e-9iy/X) ^ yl/2g-g{y/X) ^ 

Ii(Y/X,X) = yi/2(-y2 +^2)3/4^ ~ +^2)3/4^ 



1 1 

Y^~Y 



X^. (B4) 



On the other hand, if we express the outer expansion Io{Y, X) in terms of the inner variable 
y, we have 

I^{Xy,X) = ^ 

y 

(B5) 

y^ y 

We see that the 1 term outer expansion of the 2 term inner expansion is equal to the 2 term 
inner expansion of the 1 term outer expansion, in agreement with the van Dyke matching 
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principle To obtain the uniform expansion, add the inner and outer expansions, and 
subtract the overlap: 

g-9(2/) j,l/2g-9(j/) Q-Xy IX ^ 

^) = ^1/2(^2 + 1)3/4 - (^2 + 1)3/4^ ^ ^ " ? ^ 7 ^ ^ 

Q-9{y) yi/2['g-g{y) _ 



1/1/2(^2 + 1)3/4 (^2 + 1)3/4 

To obtain the small X behavior of Fq{X), we can integrate I^^if{y,X) on y. From the 
arguments given in Sec. IV. C, we know that -fo(O) = (7r/2)^/^, so the first integral is 

/■oo g-sfo) 

io ,i/2(,2 + 1)3/4 dy = vr. (B7) 
The second integral must be evaluated numerically, with the result 

-^l (^2 + 1)3/4 ^ dy = -0.7457 X. (B8) 

The last three integrals require some care. The third integral is logarithmically divergent 
for small y, so integrate down to a cutoff a, and take a — > at a later point: 

/ rfy = -(1 -7)X + Xln(a) +Xln(X) + 0(a), (B9) 

Ja y^ 

with 7 = 0.5772 .... The last two integrals are logarithmically divergent for large y (the two 
logs will cancel), so integrate up to A and then take A oo. For the fourth integral we 
then have 



X r ^ = X[\n{A) - ln(a)]. (BIO) 

Ja y 



The fifth integral converges for small y, so we can extend the lower limit of integration to 0. 
By integrating by parts, we can extract the leading In(yl) behavior; there is one remaining 
integral which must be calculated numerically, with the result 

- ^ C (^2 ^'^1)3/4 dy = -X[HA) + ln(2) - 0.4388]. (Bll) 
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Adding Eqs. (p7| ) - (|BTl| ) together, we see that the dependence upon a and A drops out, as 
it should, and we finally obtain 

2 foo Q-Xy-g{y) 



FoiX) 



dy 



(27r)l/2 7o yV2(y2 + 1)3/4 

— 4tt7 W - IA22SX + X\nX + OiX"^] 



(B12) 



This is slightly different from Eq. (37) of Larkin and Ovchinnikov [0, which in our notation 
is 



^iLarkin 




{X) 



1 



(27r)i/2 
1 



TT - ln(4e-^)X + X In X + 0(X2 
TT - 0.8091X + X InX + 0(X2 



(B13) 



(27r)i/2 

The difference between these two expressions becomes significant for values of X near 1. 
For instance, numerically we find Fq{1) = 1.732, while Eq. ( pi2|) gives 1.718, a difference of 
only 0.8%; the Larkin and Ovchinnikov expression gives 2.332, an error of 35%. The reason 
for the discrepancy is not clear, since the derivation of their result appears not to have been 
published. 



APPENDIX C: SMALL e BEHAVIOR OF THE MAGNETIZATION 

In this Appendix we will determine the small-e behavior of the dimensionless magne- 
tization, m(s), given by Eq. (|4.36| ), using the method of matched asymptotic expansions. 
The derivation is quite analogous to derivation of the small X behavior of -Fo(-^) which was 
discussed in Appendix B. Call the integrand in Eq. (|4.36|) J{y,e): 

Expand for small-e, using the small-x behavior of g{x) given in Eq. (g(y)), to obtain the 
inner expansion: 

J,(y, 6) = \5/, ^ [l + (y/vr) He/ey)e + 0{e')\ . (C2) 



31 



This expansion breaks down when = 0(1). To derive the outer expansion, define an outer 
variable Y = ey, rewrite J{y, e) in terms of Y, and expand to lowest order in e: 

e-9{Y) 



Jo{Y,e) 



(C3) 



To match the two expressions, write the inner expansion Ji{y, e) in terms of the outer variable 
Y and expand for small e: 



MY/e, e) = [1 ^ ^y/tc) ln(e/F)] e'/\ 



(C4) 



Next, take the outer expansion Jo{Y, e) and write it in terms of the inner variable y and 
expand for small e: 



Jo{ey,e)=y-'/' [1 + (y/vr) ln(e/e?/)e] . 



(C5) 



Again, we see that the 2 term inner expansion of the 1 term outer expansion is equal to the 
1 term outer expansion of the 2 term inner expansion [Q. To obtain the uniform expansion 
(i.e., an expansion valid for arbitrary y and small-e), add the inner and outer expansions, 
and subtract the overlap: 



{y + l)e-'y + y-l 



,5/2 



+ (I/tt) ln(e/ey) 



{y + l)e-^y - 1 



,3/2 



1/3/2(1 + e2|/2)3/4 



+ y 



-3/2 



(1 + e2?/2)3/4 



- 1 



(C6) 



We must now integrate Jumf{y, e) over y. For the first term we have 



°° (y + l)e~^y + y - 1 , 2(27r)i/2 
^H7^ ^^ = ^^' 



For the second term we have two integrals: 
ln(e/e) {y + l)e-'^y - 1 



(C7) 



TT Jo 



,3/2 



dy 



(27r)V' 



r^eln(e/e) 



(C8) 



e (^ + i)e-2^-l 
137^ 



TT JO 



2^ 1/2 



TT 



-4 + -ln2 + -7 

2 2 ' 



(C9) 
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The fourth and fifth integrals are both 0(e^/^), as can be seen by rescahng the integration 
variable. The fourth integral was performed numerically, with the result 

roo p-g(x) _ \ 



.1/2 



/ 



(CIO) 



x3/2(l + x2)3/4 

where the factor of 2 was determined to an accuracy of 1 part in 10^. The last integral is 

1 



Jo 



- 1 



(Cll) 



;i + x2)3/4 

We see that for all purposes the last two integrals sum to zero, although this has not been 
proven analytically. Collecting together the other terms, we finally have 



ooe-ff(^2/)[(y + l)e-22/ + y-l] 2(27r)V2 



9 

An 



In 



'8g7-5/3~ 



(C12) 



^5/2 [1 + ^2^2] 3/4 3 

A numerical evaluation of the integral for e = 1 gives 0.49577; the expansion at e = 1 gives 
0.48624, an error of about 2%. We see that the expansion is quite accurate even for relatively 
large values of e. 
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FIGURES 

FIG. 1. Illustration of the film geometry considered in this paper. The film has a width 2a in 
the x-direction, and a thickness d. The applied field Ha is in the z-direction. 

FIG. 2. Schematic diagram of the outer, inner, and matching regions used in the small-e asymp- 
totic analysis. 

FIG. 3. The inner solution F(){X, s), calculated from Eq. ( [4.30| ) (solid line). Also shown are the 
asymptotic expansions for small X (dotted line) and large X (dashed line). 



FIG. 4. Uniform approximation fumi{x,s) to the current density for e = 0.1, from Eq. ( [4.35 ) 



(solid line). Since the current density is odd in x only the current density over half the strip is 
shown. For comparison the ideal screening case (the outer solution /o(x) = — x/(l — x'^)^^'^ ) is also 
shown (dashed line). At the left edge we have fumf{x = —1) = 3.948. 

FIG. 5. Perpendicular component of the magnetic field Hz{x) within the strip for e = 0.1, 
from Eq. (^]^). The field is even in x so only half the strip is shown. Note the weak logarithmic 
singularity at x = —1. 

FIG. 6. Scaling function J^i{u) for the current density near an edge after the magnetic field is 
suddenly switched on, from Eq. ( |6.6| ). There is a maximum at u = 0.386, with J-i^ma.x = 1.1078. 



FIG. 7. Scaling functions for the response to an ac magnetic field, from Eq. (6.10); Re.7-2 is 
the solid line and lmJ-'2 is the dotted line. The real part has a maximum at n = 0.232, with 
xaax — 1.0787; the imaginary part changes sign at it — 2.43. The integral which defines the 
scaling function converges quite slowly, resulting is some numerical inaccuracies which are reflected 
in the small amplitude oscillations in the plots. 
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TABLES 

TABLE I. Summary of the primary results. The small parameter for the asymptotic expansion 
is e = 2Acg /a, where Agg is the effective penetration depth of the magnetic field and a is either half 
the width of a strip or the radius of a disk. The strip thickness is d, D = l/4Tra{0) is the diffusion 
constant for the magnetic field, and r = ad/2'iTD is a relaxation time. For a strip, ci = 0.194, 

C2 = 8.45; for a disk, ci = 0.286, C2 = 4.23. 

Outer solution in center (conformal mapping) foix) 
Inner solution near the edges (Wiener-Hopf method) Fq{X) 
Uniform solution for current density 
Uniform solution for magnetic field 
Time-dependent magnetization 
Ac magnetization 

Scaling of current at edge after jump in the field 

Velocity of current propagation after jump in the field v = 0.772D/d 
Scaling of ac current at edge uj^/'^J^2{adxoj /2D) 

Fundamental relaxation time for a strip tq = {8/10Tr)ad/D = 0.255ad/D 



/unif (-^j 'S) 

H,{x,z = 0) 

m{t) = 1 - ci(f/r) ln(1.526c2T/t) + 0{t^) 

m{uj) = 1 — ci \n{c2iijJT) / {iujT) + 0(u;~^) 
t-^l'^J^i{adx/2Dt) 
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